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Abstract 

We study a general class of quadratic BSDEs with terminal value in for p > 1. We 
prove a L^-type estimate, existence, comparison theorem, uniqueness and stability result. 

In order to construct a solution, we use a combination of the localization procedure de¬ 
veloped by Briand and Hu [6] and the monotone stability result. We point out that our 
existence result relies on rather weak assumptions on the generators. By additionally as¬ 
suming monotonicity and convexity, we deduce the comparison theorem, uniqueness and 
stability result via 0-technique developed by Briand and Hu [7]. The Krylov estimate is 
used throughout to derive some auxiliary results. We point out that our results extend 
(sometimes partially) [19], [4], [3], [8], [6], [7], [1], etc. Finally, we derive a probabilistic 
representation for the viscosity solution to the associated quadratic PDEs. 

Keywords: quadratic BSDEs, Krylov estimate, convexity, FBSDEs, quadratic PDEs 

1 Introduction 

In this paper, we are concerned with R-valued backward stochastic differential equations 
(BSDEs) 

= ^ + l F{s,Y„Z,)ds- ZsdWs, ( 1 ) 

where the generator F is continuous and satisfies P-a.s. for all (t, y, z) G [0, T] x R x 

sgn(y)F(t, y, z) < at + /5|y| -f 7^1 + /(ll/DNP, 

\F{t,y,z)\ <at + y}{\y\)+-t\z\ + f{\y\)\z\^, (2) 

for an R.+-valued progressively measurable process a, /3 G M, 7 > 0, a function /(| • j) : 

K —>■ R+ which is integrable and bounded on any compact subset of M, and a continuous 
nondecreasing function ip : R+ —>• K+. A solution to (1) is a process {Y, Z) adapted to the 
filtration generated by the Brownian motion W such that (1) holds P-a.s. for all t G [O,^]. 

We emphasize that, unlike the quadratic BSDEs studied by Briand and Hu [6], [7], the 
quadratic growth in our study takes the form /(jyDI^p. Moreover, we assume that the 
terminal value ^ and agds belong to for a certain p > 1. 
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Let us recall that, quadratic BSDEs are first studied by Kobylanski [15], where existence, 
uniqueness, comparison theorem and monotone stability for bounded solutions are obtained. 
Proving the existence of a solution consists in constructing a monotone sequence of bounded 
solutions of better-known BSDEs and then passing the limit. The underlying machinery 
of this procedure is called the monotone stability of quadratic BSDEs. Later, Briand and 
Hu [6], [7] extend the existence result by assuming that the terminal value has exponential 
moments integrability. Recently, Bahlali et al [1] constructs a solution to quadratic BSDEs 
with its terminal value in and a generator satisfying 

\F{t,y,z)\ <a + l3\y\ -f 7 |z| -f/(|?/|)|zp, 

for some a,f3,"f > 0. However, as to the uniqueness of a solution, only purely quadratic 
BSDEs are studied. 

There are two lines of studies on the uniqueness of a solution to quadratic BSDEs. 
When the terminal value is bounded, one crucial feature is that ZgdWg is a BMO mar¬ 
tingale. This observation, combined with a local Lipschitz condition, can be used to prove 
a uniqueness result; see, e.g., [13], [17[, [18], [5]. However, ZgdWs is in general not a 
BMO martingale if the terminal value is unbounded. Nevertheless one can also obtain a 
uniqueness result, by relying on a convexity condition which proves to be convenient to 
treat the quadratic generators; see [7], [14], [10], etc. 

The first contribution of this paper is to study an existence result given (2) and a 
terminal value in for a certain p > 1. We first briefly present the motivations to assume 
(2). Among the literature on non-quadratic BSDEs, assumptions of this type are quite 
convenient to obtain the a priori estimates; see, e.g., [4], [3], [8]. It turns out that the 
existence and monotone stability of bounded solutions can also be adapted to quadratic 
BSDEs with a growth of this type. This is stated in Briand and Hu [7], which assumes that 

sgn(y)F(t, y, z) < at + P\y\ + r]\z\'^, 

\F{t,y,z)\ < at + ipi\y\) + r]\z\^. 

The proof is merely a slight modification of Kobylanski [15]. In parallel with these works, 
we prove an existence result under (2). In the first step, we derive a L^-type estimate for 
quadratic BSDEs, by adapting the method developed by Briand et al [3]. To construct a 
solution, we use a combination of the localization procedure developed by Briand and Hu 
[6] and the monotone stability result. 

Another contribution is to address the question of uniqueness. In the spirit of Briand 
and Hu [7], we prove comparison theorem, uniqueness and a stability result via 0-technique 
under a monotonicity and convexity assumption. It turns out that, our results of existence 
and uniqueness, not simply provide a broader perspective in quadratic BSDEs, but also, by 
setting /(I ■ I) = 0, (partially) generalize [19], [4], [3], [8], etc. Hence our approach can be 
seen as unified to the study of both quadratic BSDEs and non-quadratic BSDEs. Einally, 
as an application, we prove a probabilistic representation for the viscosity solution of the 
quadratic PDEs associated with the BSDEs of our study. 

This paper is organized as follows. In Section 2, we introduce some functions used to 
treat the quadratic generator in (2). In Section 3, we prove the ltd-Krylov formula and a 
generalized ltd formula for y i—>■ \y\^{p > 1). The former one is used to treat discontinuous 
quadratic generators or discontinuous quadratic growth, and the later one is used to deduce 
the a priori estimates. Section 4 reviews purely quadratic BSDEs and studies their natural 
extensions, based on Bahlali et al [2]. Section 5 concerns existence, comparison theorem, 
uniqueness, etc. Finally, in Section 6, we derive the nonlinear Feynman-Kac formula in our 
framework. 
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Let us close this section by introducing all required notations. We fix the time horizon 
0 < T < +00 and a d-dimensional Brownian motion {Wt)o<t<T defined on a complete 
probability space (O, P). {J^t)o<t<T is the filtration generated by W and augmented by 

P-null sets of -P. Any measurability will refer to this filtration. In particular, Prog denotes 
the progressive cr-algebra on x [0,T]. 

As mentioned before, we only deal with M-valued BSDEs of type (1). We call the 
Prog (g)S(IR.) (8>S(M'^)-measurable random function F : 17 x [0, T] x K. x ^ K. the generator 
and the J^’-measurable random variable ^ the terminal value. The conditions imposed on 
the generator are called the structure conditions. For notational convenience, we sometimes 
write {F,^) instead of (1) to denote the BSDE with generator F and terminal value 

f^ZgdWg, sometimes denoted by Z ■ W, refers to the vector stochastic integral; see, 
e.g., Shiryaev and Cherny [22]. We call a process (Y, Z) valued in K x a solution of (1), 
if y is a continuous adapted process and Z is a Prog-measurable process such that P-a.s. 

\Zs\'^ds < - 1-00 and \F{s,Ys, Zs)\ds < -l-oo, and (1) holds P-a.s. for any t G [0,T]. 
The first inequality above ensures that Z is integrable with respect to W in the sense of 
vector stochastic integration. As a result, Z ■ W is a continuous local martingale. 

As will be seen later, the BSDEs (1) satisfying (2) is solvable if /(| • |) belongs to I, the 
set of integrable functions from K to M which are bounded on any compact subset of K. 

For any random variable or process Y, we say Y has some property if this is true 

except on a P-null subset of 17. Hence we omit “P-a.s.” in situations without ambiguity. 

Define sgn(a:) := For any cadlag adapted process Y, set Yg^t '■= Y* — Yg and 

Y* := sup(g[Q y] |Yt|. For any M-valued Prog-measurable process Ft, set |iL|s,t := Hudu 
and \H\t := jLf|o,t- T stands for the set of stopping times valued in [0, T] and S denotes the 
space of continuous adapted processes. For any local martingale M, we call {cTn}neN+ C T 
a localizing seguence if cr„ increases stationarily to T as n goes to +oo and is a 

martingale for any n G N+. For later use, we specify the following spaces under P. 

• S°°: the set of bounded processes in <S; 

• S'P{p > 1): the set of Y G 5 with Y* G L^; 

• T>-. the set of Y G <S such that {Y,-|r G T} is uniformly integrable; 

• M.\ the space of valued Prog-measurable processes Z such that P-a.s. jJ" \Zg\‘^ds < 
-l-oo; for any Z € A4, Z ■ W is a continuous local martingale; 

• M.^{p > 0): the set of Z G Ad with 


II^IUp :=E 



< -l-oo; 


• C^(R): the space of p times continuously differentiable functions from R to R; 

• FVj the Sobolev space of measurable maps u : R —>■ R such that both u and its 

generalized derivatives u',u” belong to Lj^^^(R). 


The above spaces are Banach (respectively complete) under suitable norms (respectively 
metrics); we will not present these facts in more detail since they are not involved in our 
study. We call (Y, Z) a solution of (1) if (Y, Z) belongs to S^xAi^. This definition simply 
comes from the fact that the existence holds if |^| + agds belongs to Analogously 
to most papers on R-valued quadratic BSDEs, our existence result essentially relies on the 
monotone stability result of quadratic BSDEs; see, e.g., Kobylanski [15] or Briand and Hu 

m. 
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2 Functions of Class X 


In this section, we introduce the basic ingredients used to treat the quadratic generator in 
(2). We recall that X is the set of integrable functions from R to K which are bounded on 
any compact subset of R.. 

u-^ Transform. For any / € X, define : R —>■ R and by 
u^{x):=J exp(^2y f{u)du^dy, 

:= exp (2 [ \f(u)\du\. 

^ J —00 ' 

Obviously, 1 < < + 00 . Moreover, the following properties hold by simple computations. 

Here we set u := for notational convenience. 

(i) u G C^(M)nyV^ and u"(x) = 2f(x)u'(x) a.e.; if / is continuous, then u G C^(R); 

(ii) u is strictly increasing and bijective from R to R; 

(iii) u~^ G C^(R) n Wf ;o^(R); if / is continuous, then u~^ G C^(M); 

(iv) ^ < 1 ^( 2 :)I < M\x\ and < u'{x) < M. 

Transform. For any / S X, define ^ R"*" by 

/■kl / \ 

{x) ■.= J u^~^\y)expy2 J f(u)dujdy. 

Set V := . Simple computations give 

(i) V G C^(R) n Wf ;Qg(R) and v”{x) — 2/(|a;|)|u'(a:)| = 1 a.e.; if / is continuous, then 
ueC2(R); 

(ii) v{x) > 0, sgn('(;'(a;)) = sgn(a:) and u"(0) = 1; 

(iii) 2 l^ < v{x) < and ^ < \v'{x)\ < M^\x\. 

In the sequel of our study, and exclusively stand for the above transforms associated 
with f G X. Hence in situations without ambiguity, we denote by u,v,M, 

respectively. 


3 Krylov Estimate and the Ito-Krylov Formula 


The first auxiliary result is the Krylov estimate. Later, it is used to prove an Ito’s-type 
formula for functions in C^(R) nyV^;Q^(R). This helps to deal with (possibly discontinuous) 
quadratic generators. As the second application, we derive a generalized ltd formula for 
y \y\^{P ^ 1) which is not smooth enough for 1 < p < 2. This is a basic tool to study 
hP{p> 1) solutions. 

To allow the existence of a local time in particular situations, we study equations of 
type 





F{s,Ys,Zs)ds 



ZsdWs 


(3) 


where C is a continuous adapted process of finite variation. We denote its total variation 
process by V.(C). Likewise, sometimes we denote (3) by (X, C, ^). The solution of (3) is 
defined analogously to that of (I). 

Now we prove the Krylov estimate for (3). A more complicated version not needed for 
our study can be found in Bahlali et al [1]. 
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Lemma 1 (Krylov Estimate) Consider (3). For any measurable function '0 : R —>■ R+, 


E 



f^{Ys)\Zs\^ds 




where is a stopping time defined by 


:=inf{t>0: \Yt\ + VtiC) + |E(s, F,, Z«)|ds > m} A T. 


(4) 


Proof. Without loss of generality we assume ||'0llLi([-m m]) ^ +oo- For each n £ N+, set 


Tm,n ■= Tm A iuf 


> 0 : J \Zs\'^ds > n|. 


Let a £ [—m,m]. By Tanaka’s formula, 


- a)“ = (4^0 - a)” - 


I 


/o 


{n<a}ars 


dYs + Ar^JY) 


pt/\Tm,n pt/\Tm,n 

= {^0 ~ + / ^{Ys<a}Fi^^^ST 2s)ds / ^{Ys<a}(^Cs 


/O 


/o 


ctATrr 


1 




(5) 


where L°^[Y) is the local time of F at a. To estimate the local time, we put it on the left-hand 
side and the rest terms on the right-hand side. Since x i—>■ (x — a)“ is Lipschitz-continuous, 
we deduce from the definition of that 


(Fo - a)- - - a)- < |Fo - Yt^r^J < 2m. 

Meanwhile, the definition of Tm also implies that the sum of the ds-integral and dC-integral 
is bounded by m. Hence, we have 


E[L“^,^ JF)] <6m. 

By Fatou’s lemma applied to the sequence indexed by n, 

sup E[L“^.,.^(F)] < 6m. 

a€[—m,m] 


We then use time occupation formula for continuous semimartingales (see Chapter VI., 
Revuz and Yor [21]) and the above inequality to obtain 


E 


nT ATjt 


f;{Ys)\Zs\'^ds 


'-JO 


= E / f^ix)L^r^,jY)dx 

— m 

f^{x)E[Lif^,jY)]dx 

< 6TO||V’llLi([-m,m])- 


□ 
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As an immediate consequence of Lemma 1, we have P-a.s. 


T 

= 0 , ( 6 ) 

for any A C K with null Lebesgue measure. This will be used later several times. 

Given Lemma 1, we turn to the main results of this section. The following generalized 
ltd formula is proved in Bahlali et al [1]. 

Theorem 2 (Ito-Krylov Formula) If {Y, Z) is a solution of (3), then for any u € C^(R)n 
we have P-a.s. for all t G [0,T], 

u(Yt) = u(Yo) + [ u'iY,)dY, + 11 u"iY,)\Zs\^ds. (7) 

Jo ^ Jo 


Proof. We use defined in Lemma 1 (Krylov estimate). Note that Tm increases station- 
arily to T as m goes to -|-oo. It is therefore sufficient to prove the equality for u{Yt/\T^). 
To this end we use an approximation procedure. We consider m such that P-a.s. m> ITqI- 
Let Un be a sequence of functions in (R) satisfying 

(i) Un converges uniformly to u on [—m,m]; 

(ii) u'^ converges uniformly to u' on [—m, m]; 

(hi) u" converges in L^([—m,m]) to u". 

By Ito’s formula, 

ptATm 1 ptATrn 

Un{Yt/\r„,)=Un{Yo)+ Un{Ys)dYs + - u'f{Ys)\Zs\'^ds. 

Jo ^ Jo 

Due to (i) and iLtAr™! < rn, Un{Y.^T-m) converges to ufY.^rm) P-a.s. uniformly on [0,T] as 
n goes to -|-c»; the second term converges in probability to 

u{Ys)dYs 

by (ii) and dominated convergence for stochastic integrals; the last term converges in prob¬ 
ability to 

- / u"{Ys)\Zs\^ds 

^ Jo 

due to (iii) and Lemma 1. Indeed, Lemma 1 implies 



E 




< - u"m)\Z^\^ds < 6m||< - 


Hence collecting these convergence results gives (7). By the continuity of both sides of (7), 
the quality also holds P-a.s. for all t G [0,T]. 


□ 

To study LP(p > 1) solutions we now prove an Ito’s-type formula for y i—;• \y\^{p > 1) 
which is not smooth enough for 1 < p < 2. The proof for multidimensional ltd processes can 
be found, e.g., in Briand et al [3]. In contrast to their approach, we give a novel and simpler 
proof for BSDE framework but point out that it can be also extended to ltd processes. 
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Lemma 3 Let p > 1. If {Y, Z) is a solution of (3), then we have P-a.s. for all t G [0, T], 

\Yt\P + l\[Y.^o}\Ysr^\Zs\^ds 

= \e-p|\gn{Ys)\Ys^^dYs-l^p=^ dL°{Y), (8) 

where L^{Y) is the local time ofY at 0. 

Proof, (i). p = 1. This is immediate from Tanaka’s formula. 

(ii). p > 2. y ^ \yY’ G C^(R). Hence this is immediate from Ito’s formula. 

(hi), p = 2. p I—>■ \yY' G C^(]R). Due to (6), /p \Ys\^~^\Zs\^da is indistinguishable from 
J^I[Y,^oy\Ys\P~^\Zs\^ds. By taking this fact into account, this equality is thus immediate 
from Ito’s formula. 

(iv). 1 < p < 2. We use an approximation argument. Define 

ue{y) ■= + e^)"- 

Hence for any e > 0, we have uf G C^(K.). By Ito’s formula, 

<(rt) = um-P[ Y,ur\Y,)dY, - {puP-\Y,)+p{p-2)\Y,\^ul-\Y,))\Z,\Hs. 

(9) 


Now we send e to 0. Ue(p) —>■ \y\ pointwise implies u^fYi)^ —)> \Yt\P and —S' 

1^1^ pointwise on D. Secondly, yu^~‘^{y) —^ sgn(p)|p|^“^ pointwise implies by dominated 
convergence for stochastic integrals that 


i;sgn(ysX '^{Ys)dYs- 


\Y,\P-~^dYs in probability. 


To prove that the ds-integral in (9) also converges, we split it into two parts and argue their 
convergence respectively. Note that 

puP-\Y,)+pip-2)\Y,\^uP-%Y,)=pe^uP-\Ys)+p{p-l)\Ys\^uP-HY,). ( 10 ) 

For the second term on the right-hand side of (10), we have 


iniV"(>^.) = im/o}inr 


IWI 


eiYs) 


4—p 


Since monotone convergence gives 


f \Ys\'^uP ^{Ys)\Zs\‘^ds —)> / ^Ipointwise on D. 

Jt Jt 


It thus remains to prove the ds-integral concerning the first term on the right-hand side 
of (10) converges to 0. To this end, we use Lemma 1 (Krylov estimate) and the same 
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localization procedure. This gives 


E 


e^uP-'^iYs)\Zs\^ds <6me‘^ {x^ + e‘^)'^dx 


'-Jo 


<12me^ (x^+e ^)''2 dx 

Jo 

pm 

< 12 • 2“^ me^ / (x + eY~'^dx 
Jo 


^m+e 


<12-2 2 me 


vP-'^dx 


12 • 2^m 
p-3 


{e^{m + e)P-^ -eP-^), 


which, due to 1 < p, converges to 0 as e goes to 0. Hence e'^uP~'^{Ys)\Zs\^ds converges 
u.c.p to 0. Collecting all convergence results above gives (8). Finally, the continuity of each 
term in (8) implies that the equality also holds P-a.s. for all t G [0,T]. 


□ 


4 > 1) Solutions of Purely Quadratic BSDEs 


Before turning to the main results of this paper, we partially extend the existence and 
uniqueness result for purely quadratic BSDEs studied by Bahlali et al [1]. Later, we present 
their natural extensions and the motivations of our work. These BSDEs are called purely 
quadratic, since the generator takes the form F{t,y,z) = /(y)|zp. The solvability simply 
comes from the function defined in Section 2 which transforms better known BSDEs to 
(/(y)kP,Q by Ito-Krylov formula. 

Theorem 4 Let f Gl and ^ £ ]IP'(p > 1). Then there exists a unique solution of 

Yt=^ + J^ f{Ys)\Zs\^ds-ZsdWs. ( 11 ) 

Moreover, if p > 1, the solution belongs to Sp x Mp; if p = 1, the solution belongs to 
V X for any q G (0,1). 

Proof. Let u := and M := . Then u,u~^ £ C^(K) n The existence and 

uniqueness result can be seen as a one-on-one correspondence between solutions of BSDEs. 

(i). Existence. |u(a;)| < M\x\ implies u(^) £ L^. By Ito representation theorem, there 
exists a unique pair (Y, Z) which solves (0, u{f)), i.e., 

dYt = ZtdWt,YT = u{0. (12) 

We aim at proving 

(y,Z) :=(„-■(?), ,, (13) 

u'{u (Y)) 

solves (11). Ito-Krylov formula applied to Yt = u~^{Yt) yields 


dYt = 


u'{u-^{Yt)) 


-dYt- 


1/ 1 \2n" 

2 7// 


2u"(u-i(y*)) 


u'iu-^{Yt))^ u'{u-^{Yt)) 


\Zs\‘^ds. 


(14) 









To simplify (14) let us recall that u"{x) = 2f{x)u'{x) a.e. Hence (13), (14) and (6) give 

dYt = -f{Yt)\Zt\^dt + ZtdWu Yt = 


i.e., {Y,Z) solves (11). 

(ii). Uniqueness. Suppose {Y,Z) and {Y',Z') are solutions of (11). By Ito-Krylov 
formula applied to u{Y) and u(Y'), we deduce that (u(Y),u'(Y)Z) and {u{Y'),u'(Y')Z') 
solve (0,u(^)). But from (i) it is known that they coincide. Transforming u{Y) and u(Y') 
via the bijective function u~^ yields the uniqueness result. 

(hi). We prove the estimate for the unique solution {Y, Z). For p > I, Doob’s hP^p > 1) 
maximal inequality used to (12) implies {Y, Z) G x MP. Hence (Y, Z) G x due 
to |u'(a:)| > ^ and |u“^(a:)| < M|a:|. For p = 1, Y gV since it is a martingale on [0,T]. 
In view of the above properties of u we have Y gT). The estimate for Z is immediate from 
Lemma 6.1, Briand et al [3] which is a version of L^(0 < p < 1) maximal inequality for 
martingales. 


□ 

Remark. If ^ is a general J^y-measurable random variable, Dudley representation theorem 
(see Dudley [11]) implies that there still exists a solution of (12) and hence a solution of 
(11). However, the solution in general is not unique. 

The proof of Theorem 4 indicates that / being bounded on compact subsets of K is not 
needed for the existence and uniqueness result of purely quadratic BSDEs. 

Proposition 5 (Comparison) Let f,g G I, G LP^p > 1) and (Y,Z), (Y',Z') be the 
unique solutions of {f (y)\z\^, , {g(y)\z\^respectively. If f < g a.e. andP-a.s. f < ff, 

then P-a.s. Y. < Y!. 

Proof. Again we transform so as to compare better known BSDEs. Set u := . For any 

T G T, Ito-Krylov formula yields 

u{Yf^,) = u{Y;) + f (u'{Yf)g{Yf)\Z'f -\u''{Yf)\Z'f)ds - f u'lY,)Z'JWs. 

JtAT ^ Z / JtAr 

= uiY;)+ r u'iYf){g{Y:)-f{Y:))\Z'fds- f u\Y,)ZidWs 

JtAr J tAr 

> u{Yf) - [ ufYs)Z'JWs, 

J tAr 

where the last two lines are due to u"{x) = 2f{x)u'{x) a.e., g > f a.e. and (6). In the next 
step, we want to eliminate the local martingale part by a localization procedure. Note that 
u'{Ys)ZgdWs is a local martingale on [t,T]. Set {T„}„gpj+ to be its localizing sequence 
on [t,T]. Replacing r by in the above inequality thus gives P-a.s. 

u{Yl)>E[u{Yf^,J\A]. 

This implies that, for any A G we have 

E[u{Yf)lA] >E[u(y,;,JI^]. 

Since u{Y') G V, we can use Vitali convergence theorem to obtain 

E[u(y/)I^] >E[u(OIa] =E[E[u(C')|-^*]Ia]- 
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Note that this inequality holds for any A G jFf Hence, by choosing A = {uiYD < 
E[u(^')| J't]}, we obtain u{Yl) > E[u(^')|j^t]. Since ^ and u is increasing, we fur¬ 
ther have uiXl) — Let us recall that, by Theorem 4, {u{Y),u'(Y)Z) is the 

unique solution of (0,u(^)). Hence, u{Yl) > u{Yt). Transforming both sides via the bijec- 
tive increasing function u~^ yields P-a.s. Yt < Y/. By the continuity of Y and Y' we have 
P-a.s. Y<Y!. 


□ 


Remark. In Proposition 5, we rely on the fact that P-a.s. 

I - /(n')w'(H;)) \z'fds = o, (is) 

even though u"{x) = 2f{x)u'{x) only holds almost everywhere on K. Here we prove it. Let 
A be the subset of K on which u"{x) = 2f{x)u'{x) fails. Hence, 


IfV/GRVA} 




iz'pds = 0. 


Meanwhile, by ( 6 ), we have P-a.s. 


I 


{^' 6^1 




iZ'l^ds = 0. 


Hence, (15) holds P-a.s. This fact also applies to Theorem 4 and all results in the sequel of 
our study. 


Theorem 4 and Proposition 5 are based on a one-on-one correspondence between so¬ 
lutions (respectively the unique solution) of BSDEs. Hence it is natural to generalize as 
follows. Set f G I, u := , F{t, y, z) := G{t, y, z) + f{y)\z\'^ and 

F(t,y,z) :=u'iu-\y))G{t,u-\y), ^ ). (16) 

u'[u Ay)) 

If G ensures the existence of a solution of {F, u{^)), we can transform it via u~^ to a solution 
of (F,^). An example is that G is of continuous linear growth in (y, z) where the existence 
of a maximal (respectively minimal) solution of (A, u(^)) can be proved in the spirit of 
Lepeltier and San Martin [16]. 

When the generator is continuous in (y, z), a more general situation is linear-quadratic 
growth, i.e.. 


\H{t,y,z)\ < a + Ply] + Az] + f{\y\)\zf :=F{t,y,z), (17) 

for some a, /3 ,7 > 0. The existence result then consists of viewing the maximal (respectively 
minimal) solution of (F,y) (respectively {—F,—^~)) as the a priori bounds for solutions 
of (iL, ^), and using a combination of a localization procedure and the monotone stability 
result developed by Briand and Hu [ 6 ], [7]. For details the reader shall refer to Bahlali et 
al [ 1 ]. 

However, either an additive structure in (16) or a linear-quadratic growth (17) is too 
restrictive and uniqueness is not available in general. Considering this limitation, we devote 
Section 5 to the solvability under milder structure conditions. 
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5 > 1) Solutions of Quadratic BSDEs 

With the preparatory work in Section 2, 3, 4, we study > 1) solutions of quadratic 
BSDEs under general assumptions. We deal with the quadratic generators in the spirit of 
Bahlali et al [1], derive the a priori estimates in the spirit of Briand et al [3] and prove the 
existence and uniqueness result in the spirit of Briand et al [6], [7], [8]. This section can 
also be seen as a generalization of these works. The following assumptions on {F, ensure 
the estimates and an existence result. 

Assumption (A.l) Let p > 1. There exist /3 G R, 7 > 0, an R^-valued Prog-measurable 
process a, /(| ■ |) G I and a continuous nondecreasing function (p : R.+ —>• R.+ with (/3(0) = 0 
such that 1^1 -I- \a\T G and P-a.s. 

(i) for any t G [0, T], (y, z) 1 — F{t, y, z) is continuous; 

(ii) F is “monotonic” at y = 0, i.e., for any (t, y, z) G [0, T] x R x R'^, 

sgn(y)F(t, y, z) < at + P\y\ + -i\z\ + f(]y\)\z\^- 

(iii) for any (t, y, z) G [0, T] x R x R"^, 

\F{t,y,z)\ < at + Lp{\y\)+-i\z\ + f(\y\)\z\^. 

It is worth noticing that, given (A.l)(iii) and /(| • |) = 0, (A.l)(ii) is a consequence of F 
being monotonic at y = 0. Indeed, 

sgn(y - 0){F{t,y,z) - F(t,0,z)) < /3\y\ 


implies 


sgn(y)F(t, y, z) < F{t, 0, z) + /3|y| 
< Oit + (3\y\ +j\z\. 


This explains why we keep saying that F is monotonic at y = 0, even though y also appears 
in the quadratic term. Secondly, our results don’t rely on the specific choice of ip. Hence 
the growth condition in y can be arbitrary as long as (A.l)(i)(ii) hold. Assumptions of this 
type for different settings can also be found in, e.g., [4], [3], [7]. Finally, /(| • |) can be 
discontinuous; /(| ■ |) being R^-valued appears more naturally in the growth condition. 


Lemma 6 (A Priori Estimate (i)) Let p > 1 and (A.l) hold for (E, Q. If (Y, Z) G 
SP X M is a solution of (E, C); then 


E 


( l^\Z.,\^dsy] f{\Ys\)\Zs\^dsy] <c(E[(W)^’ + |a|?,]) 


Tj > 


where c is a constant only depending on 

Proof. Set V := and M := For any t Ito-Krylov formula yields 


;(yo) = viYr) + f v{Ys)F{s, Ys,Zs)ds 

Jo 

-i/ v"{Ys)\Z,\^ds- [ v'{Ys)Z,dW,. 
^ Jo Jo 


(18) 
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(19) 


Due to sgn(w'(a;)) = sgn(a:) and (A.l)(ii), we have 

v'{Y,)F{s,Y,,Z,) < \v'{Y,)\{at+m\Yl\Zs\+f{\Y,\)\Z,\^). 

Recall that v"{x) — 2f(\x\)\v'{x)\ = 1 a.e. Hence (18) and (19) give 

i / \Zs\^ds <v{Yr) + [ \v'{Ys)\{as+l3\Ys\+^\Zs\)ds- f v'(Ys)ZsdWs. 

^ Jo Jo Jo 

Moreover, since v(x) < and |u'(a;)| < M^\x\, this inequality gives 

f \Zsfds<ci{Y*f + cif \Y,\{as + \Ys\ + \Zs\)ds-2 f v'{Ys)Z,dWs, 

Jo Jo Jo 

where ci := 2M^(1 V ,5 V 7 ). Note that in (20), 

i(y *)2 + i|a||, 

ci^ \Ys\\Zs\ds<^clT-{Y*f + ^J^ \Z,fds. 

Hence (20) yields 

f |Z,| 2 ds< (3ci+c?T)(r*)2+ci|a||-4 / v'iYs)Z,dW,. 

Jo Jo 

This estimate implies that for any p > 1, 


( 20 ) 


E 




< C 2 E 


(y*) 


*'\P _L U|P 


v'iYs)ZsdWs 


( 21 ) 


where C 2 := 3= ((3ci + cfT) V 4) ^. Define for each n G N+, := inf {f > 0 : f* \Zs\^ds > 

n} A T. We then replace t by r„ and use Davis-Burkholder-Gundy inequality to obtain 


C 2 E 


(^j\'{Y,)Z,dW,y^ < C2c(p)MPe[(^ ” \Ys\^\Z,\^dsy 


< \clc{pfM^^¥.[{Y*Y] + ifi 

< + 00 . 




We explain that in this inequality, c(p) comes from Davis-Burkholder-Gundy inequality and 
only depends on p. With this estimate, we come back to (21). Transferring the quadratic 
term to the left-hand side of (21) and using Fatou’s lemma, we obtain 


E 



£ 

2 


<c(E[(F*)P + |a|^]), 


where c := c^YpYM'^p -|- 2 c 2 . 

Toestimate /(iWDlZsI^ds we use u := This helps to transfer f{\Ys\)\Zs\'^ds 

to the left-hand side so that standard estimates can be used. The proof is omitted since it 
is not relevant to our study. 


□ 
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We continue our study by sharpening Lemma 6 for p > 1. We follow Proposition 3.2, 
Briand et al [3] and extend it to quadratic BSDEs. As an important byproduct, we obtain 
the a priori bound for solutions which is crucial to the construction of a solution. 


Lemma 7 (A Priori Estimate (ii)) Let p > 1 and (A.l) hold for (F,^). If {Y,Z) € 
X M is a solution to (E, ^), then 

E[(W)?’] +e[(^I^ \Zs\^dsY] +e[(^I^ fi\Ys\)\Zsfdsy] <c(E[|^|^’ + |a|?,]). 

In particular, 


E 


sup \Y,\P Ft <cE[|^|P + |a|P^|j-*]. 

LsGft.T] 


In both cases, c is a constant only depending on T, /3, 7 ,p. 

Proof. Let u := and M := and denote wdYtl), u'{\Yt\), u"{\Yt\) by ut, u't, u'f, 

respectively. By Tanaka’s formula applied to \Yt\ and Ito-Krylov formula applied to ut, 

Ut = UT + J^ sgn{Ys)u',,F{s,Ys,Zs)ds-^ l[Y,^o}u's\Zs\'^ds 

- [ sgn{Ys)u,,ZsdWs - [ uJL°{Y), 

Jt Jt 

where L^iY) is the local time of Y at 0. Lemma 3 applied to ut then gives 
= \ut\^ +p sgn(us)|usr"^(sgn(WXE(s, 

-p f sgn{us)\us\^~^u',,dL°{Y) - p f sgn(Ms) sgn(W)|usr“^’^s^sclW4- 

Jt Jt 

To simplify this equality, we recall that sgn(Ms) = = I{y, 5 ^o} und u"{x) = 2f{x)u'{x) 

a.e. Hence 

\ut\^ + \Y.^O}\UsrW\Zs?ds 

< \uTf’+p I{Y,^o}\us\^~^u',,{as + l3\Ys\ + -f\Zs\)ds 

-pj^ sgn(W)|wsr"^w's^sdbTs- 

Let {c„}„gN+ be constants to be determined. Since < u(|a;|) < M\x\ and -F < M'(|a;|) < 
M, this inequality yields 

< + mYs\+l\Zs\)ds 

-pj^ sgn{Ys)\us\P~^u',,ZsdWs, ( 22 ) 
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where ci := ^ Observe that in (22), 

MPj^y^^o}\ysr^\Zs\ < 

We then use this inequality to (22). Set C 2 := Mp V (Mp|/3| + ^ 2 ci )’ 

^ ■= ^2(1^1'’ + inr'K + \Ys\)ds), 

and N to be the local martingale part of (22). Hence (22) gives 

|Ftr + y^ IiY^^o}\Ysr^\Zsfds<X-NT + Nt. (23) 

We claim that TV is a martingale. Let c(l) be the constant in Davis-Burkholder-Gundy 
inequality for p = 1. We have 


E[N*] < c(l)E[(iV)' 


< c(l)MPE 

^ c(l)MP 
P 

< +00, 


\Ys\^P-^\Z,\^ds^ 

((p-i)E[(r*)p] +e[( y" \Zs\^ds 


where the last two lines come from Young’s inequality and Lemma 6 (a priori estimate (i)). 
Hence iV is a martingale. Coming back to (23), we deduce that 


E 




< —Em. 

Cl 


(24) 


Now we estimate Y via X. To this end, taking supremum over t G [0,T] and using Davis- 
Burkholder-Gundy inequality to (23) give 


E[(r*)P] <E[Y]-kc(l)E[(Y)|]. 


(25) 


The second term in (25) yields by Cauchy-Schwartz inequality that 


c(l)E[(7V)|:] <c(l)MPE (Y*)i( / 


< 2 E[(Y*)^’] 


c( 1 ) 2 m 2 p 


E 


'-^0 


l{Y.^o}\Ysr^\Zs\^ds 


Using (24) to this inequality gives the estimate of {N)^ via Y and X. With this estimate 
we come back to (25) and obtain 


E[(r*)P] < 2 ^ 


1 -k 


2c(1)2m2p 


Cl 


)e[X]. 


Set C 3 := 2 c 2(1 -I- ^ This inequality yields 


E[(Y*)P] <C 3 (E[|er] +E[y \Ys\P-^asds\+E[J^ I^Rs]) 


(26) 
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Young’s inequality used to the second term on the right-hand side of this inequality gives 


C3 r \Y,r^asds < l(Y*r + A ) , 

Jo 2 p Vcsg/ 

where q is the conjugate index of p. Set C 4 := 2^03 V (26) and the above 

inequality yield 


mY*r] <c4(E[iei' 


-E 


■- ^0 


sup 

u.G[ 0 ,s 




By Gronwall’s lemma, 


E[{Y*f] < C4exp(c4T)E[|^|P + |a|P]. 


Finally, by Lemma 6 we conclude that there exists a constant c only depending on T, M, /3, 7 , p 
such that 


E[(y*)p] -hE 


iZ.J'^ds] 


-l-E 


fi\Ys\)\Zs\^dsy] <cE[|^r + |a|?,]. 


To prove the remaining statement, we view any fixed t G [OiT] as the initial time, reset 


^ := C2(ier +inr'K + mods) 


and replace all estimates by conditional estimates. 


□ 


An immediate consequence of Lemma 7 is that 

\Yt\ < (cE[|er + |a|^|J■,])^ 

i.e., Y has an a priori bound which is a continuous supermartingale. 

With this estimate we are ready to construct a LP(p > 1) solution via inf-(sup-)convolution 
as in Briand et al [ 6 ], [7], [ 8 ]. A localization procedure where the a priori bound plays a 
crucial role is used and the monotone stability result takes the limit. 

Theorem 8 (Existence) Let p > 1 and (A.l) hold for (F,C)- Then there exists a solution 
of {F,C) inSP X AdP. 

Proof. We introduce the notations used throughout the proof. Define the process 

At := (cE[|Cr + |o|?.|At])^ 

where c is the constant defined in Lemma 7. Obviously X is continuous by Ito representation 
theorem. Moreover, for each m,n G N'*', set 

Tm ■= inf {t > 0 : |a|t -I- At > m} A T, 

(T„ := inf {t > 0 : jolt > n} A T. 
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It then follows from the continuity of |a|. and X that and (t„ increase stationarily to T 
as TO,n goes to +oo, respectively. To apply a double approximation procedure, we define 


y,z) := l{t<a^} inf {F+{t, y', z') + n\y 




“ + fell/ 


- y'\ + n\z - z'\} 
-y'\ +k\z- z'\}, 


and := A n — A /c. 

Before proceeding to the proof we give some useful facts. By Lepeltier and San Martin 
[16], is Lipschitz-continuous in {y,z); as k goes to +oo, converges decreasingly 
uniformly on compact sets to a limit denoted by as n goes to +oo, F”’°° con¬ 

verges increasingly uniformly on compact sets to F. Moreover, | |F”’^(-, 0, 0)| and 
are bounded. 

Hence, by Briand et al [3], there exists a unique solution G <5^ x of 

by comparison theorem, is increasing in n and decreasing in k. We are 
about to take the limit by the monotone stability result. 

However, | |F”’^(-, 0, 0)| and Y'^’^ are not uniformly bounded in general. To over¬ 
come this difficulty, we use Lemma 7 and work on random time interval where and 
I |F"’^(-, 0,0)11 are uniformly bounded. This is the motivation to introduce X and t^. To 
be more precise, the localization procedure is as follows. 

Note that (F”’*^, satisfies (A.l) associated with (a,/3, 7 , </?,/). Hence by Lemma 7 
(a priori estimate (ii)), 

< A*. (27) 

In view of the definition of Tm, we deduce that 

jy^vtl < <m. (28) 

Hence is uniformly bounded on [0,Tm]- Secondly, given which solves 

it is immediate that (Y:”vt, I[o.r,„]2'”’'=) solves (I[o,^^]F"’^ To make the 

monotone stability result adaptable, we use a truncation procedure. Define 

Piy) I{y< —m}^ F ^{\y\<m}y “b 

Hence from (28) (K”’.,^, I[o,t„]^"’^) meanwhile solves (I[o_T-^](t)F”’''(t, p(j/), z), Sec¬ 

ondly, we have 

ho,r^]it)F^’'^{t,p{y),z)\ < (at + (/7(|p(y)|) +7|z| + /(|p(y)|)|zp) 

<ht<Tr,,}(at + pim)+j\z\+ sup /(|p(?/)|)|zp) 

^ lyl^m ^ 

2 

<ht<r,^}(at + p{m)+ '^ + { sup f{\p{y)\) + l)\z\A, 

where sup|j,|<^ fi\piy)\) is bounded for each m due to /(| • j) G X. Moreover, the definition 
of Tm implies \a\r^ < m. Hence we can use the monotone stability result (see Koby- 
lanksi [15] or Briand and Hu [7]) to obtain (y"*.".oc^ g x which solves 

(I[Q (t)F”’°°(t, p(j/), z), inffe YJ^^). Moreover, is the P-a.s. uniform limit of 

as k goes to + 00 . These arguments hold for any m,n £ N"*". 
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Due to this convergence result we can pass the comparison property to We use 

the monotone stability result again to the sequence indexed by n to obtain Z™) € 

X which solves (I[o,t-„] p(y), 2 :), sup„ inffc Likewise, F™ is the P-a.s. 

uniform limit of as n goes to + 00 . Hence we obtain from (28) that |F(™| < XtAr^ < 

m. Therefore, (V™,Z"^) solves (I[o_t-^]F, sup„ inf^ i.e.. 


F 


tATrr 


supinfF,’:;;'=+ / Fis,Yr, ZT)ds - Z^dWs- 

^ ^ JtATm JtATm 


(29) 


We recall that the monotone stability result also implies that F"* is the A^^-limit of 
as fc,n goes to +c>o. This fact and previous convergence results give 

F™+i = F-P-a.s., 

= I{t<^„}FF dt ® dP-a.e. (30) 

Define (F F) on [0,r] by 


F := ip<.aF' + E 

m>2 

F := ht<r,}Zl + ^ Ip„,_„.^]Fr. 


'i>2 


By (30), we have = Y^^^ and Ip<,-^}Ft = Ip<,-^}FF. Hence we can rewrite (29) as 

F(s,F,F«)ds- f ” ZsdWs. 

By sending m to + 00 , we deduce that (F, F) solves {F, ^). Since (F"’^, F"’^) verifies Lemma 
7, we can use Fatou’s lemma to prove that {Y,Z) G x Ad^. 


F 


tATr, 


= sup inf FF’" + 


□ 

Theorem 8 proves the existence of a L^’(p > 1) solution under (A.l) which to our 
knowledge the most general asssumption. For example, (A.l)(ii) allows one to get rid of 
monotonicity in y which is required by, e.g., Pardoux [19] and Briand et al [4], [3], [8]. 
Meanwhile, in contrast to these works, the generator can also be quadratic by setting 
/(I • I) S I. Hence Theorem 8 provides a unified way to construct solutions of both non¬ 
quadratic and quadratic BSDEs via the monotone stability result. 

On the other hand. Theorem 8 is an extension of Bahlali et al [1] which only studies 
BSDEs with L2 integrability and linear-quadratic growth. However, in contrast to their 
work, (A.l) is not sufficient in our setting to ensure the existence of a maximal or minimal 
solution, since the double approximation procedure makes the comparison between solutions 
impossible. 

However, to prove the existence of a maximal or minimal solution is no way impossible. 
Since we have X as the a priori bound for solutions, we can convert the question of existence 
into the question of existence for quadratic BSDEs with double barriers. This problem has 
been solved by introducing the notion of generalized BSDEs; see Essaky and Hassani [12]. 

Let us turn to the uniqueness result. Motivated by Briand and Hu [7] or Da Lio and 
Ley [9] from the point of view of PDEs, we impose a convexity condition so as to use 
0-techiiiique which proves to be convenient to treat quadratic generators. We start from 
comparison theorem and then move to uniqueness and stability result. To this end, the 
following assumptions on (F, ^) are needed. 
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Assumption (A.2) Let p > 1. There exist /3i,/32 S R, 71 j 72 > 0, an R+-valued Prog- 
measurable process a, a continuous nondecreasing function ip : M+ R+ with (p(0) = 0, 
/(I • I) el and Fi,F 2 : x [0, T] x K x —>• K which are Prog 0 B(R'^)-measurable 

such that F = Fi + F 2 , |^| -I- |q!|t S and P-a.s. 

(i) for any t G [0, T], {y, z) 1 — >■ F{t, y, z) is continuous; 

(ii) Fi{t, y, z) is monotonic in y and Lipschitz-continuous in z, and p 2 (t, y, z) is monotonic 
at 2 / = 0 and of linear-quadratic growth in 2 , i.e., for anyt G [ 0 ,T],y,y' G R, z' G 

sgn{y - y'){Fi{t,y,z) - Fi{t,y',z)) < / 3 i|y - y'\, 

\Fiit,y,z) - Fi{t,y,z')\ < 7i|z - z'\, 

sgn(y)F2(t,y,z) < /Jalyl + 72^1 -h/(|y|)|zp; 

(iii) for any t G [0, T], {y, z) 1 — >■ F2{t, y, z) is convex; 

(iv) for any (t, y, z) G [0, T] x R x R^', 

\F{t,y,z)\ <at + p{\y\) + (71 +72)1^1 + /(ll/DkP- 


Intuitively, (A.2) specifies an additive structure consisting of two classes of BSDEs. 
The cases where F 2 = 0 coincide with classic existence and uniqueness results for R-valued 
BSDEs; see, e.g., Pardoux [20] or Briand et al [4], [3]. When Fi = 0, the BSDEs concern and 
generalize those studied by Bahlali et al [2]. Given convexity as an additional requirement, 
we can prove an existence and uniqueness result in the presence of both components. This 
can be seen as a general version of the additive structure discussed in Section 4 and a 
complement to the quadratic BSDEs studied by Bahlali et al [2] and Briand and Hu [7]. 

We start our proof of comparison theorem by observing that (A.2) implies (A.l). Hence 
the existence of a L^(p > 1 ) solution is ensured. 

Theorem 9 (Comparison) Let p > 1, and {Y, Z), {Y', Z') G x A4 be solutions of 
{F,f), {F'^f), respectively. Iff’-a.s. for any {t,y, z) G [0,T]xRxR‘^, F{t,y,z) <F'(f,y,z) 
and f and F verifies (A.2), then P-a.s. Y. < Y.'. 


Proof. We introduce the notations used throughout the proof. For any 6 G (0,1), define 

SFt ■.= F{tX,Zi)-F'{tX,Zf), 

SsY := y - BY', 

SY ■.= Y - Y', 


and 5gZ,5Z, etc. analogously. 6 *-technique applied to the generators yields 


F{t,YuZt)-eF'{t,Y;,Z'^) 

= {Fit, Yt,Zt) - 9F{t, Yf, Z{)) + 0{F{t, Yf, Z() - F'{t, Yf, Z')) 

= eSFt + {Fit, Yt, Zt) - BF{t, Yf, Z()) 

= 95Ft + (Fi(t, Yt, Zt) - 9F^{t, Yf, Z'f)) + {F2{t, Yt, Zt) - 9F2{t, Yf, Z'^)). 


By (A. 2 )(iii), 


F2{t,Yt,Zt) = F2{t,9Yf + (f - 0)^,9Z't + (f - 9)^) 


<9F2{t,Yf,Z't) + {l-9)F2{t, 


SgYt 6gZt 
1-9’ 1-9 


)■ 


(31) 
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Hence we have 


F 2 {t, Yu Zt) - 0F2it, y;, ZO < (1 - 0)F,it, ^). (32) 

Let u be the function defined in Section 2 associated with a function of class I to be 
determined later. Denote u{{ 6 gYt)'^),u'{{d 0 Yt)'^),u"{{ 6 gYt)'^) by ut,u[, u", respectively. It 
is then known from Section 2 that Ut > 0 and u[ > 0. For any t € T, Tanaka’s formula 
applied to (SgY)'^, Ito-Krylov formula applied to u{{SgYt)'^) and Lemma 3 give 




P{P - 1) 


I{S,Y.>O}\UsrW\S0Zs\'^ds 


< \u.\P 


+ p r I|5,y.>0}|Usr' U{F{s,Ys,Z,)- 0 F'{sX.Zi)) - \u'XoZs?) ds 

JtAr si__ £ _ 


:=A. 


-p [ ^u'JgZsdWs. 

J tAr 


(33) 


By (31), (32), (A.2)(ii) and SF < 0, we deduce that, on {SgYg > 0}, 


A, < 


f(Miyi) . I 

/,(^FX.Ys,Zs)-0FAsX,Z:) + P2{SgYs)++j2\SgZs\ + \SgZ,\^) - -u'XeZs?. 




To eliminate the quadratic term, we associate u with Ag i 


u{x) : = / exp 


H 


V ffJiiL') . 
^^du)dy 


= exp ( 2 / f{\u\)du^dy. 

Hence, on {SgYg > 0}, the above inequality gives 

< <(Fi(s, Ys,Zs) - 0 FAs, y;, Z() + XddYX + -i 2 \ 5 gZs\). 


(34) 


We are about to send 0 to 1, and to this end we give some auxiliary facts. Reset M := 
exp (2 JA f{u)du). Obviously 1 < M < + 00 . By dominated convergence, for x > 0, we 
have 


lim u(x) = Mx, 
e^i 

lim u'{x) = MI{a,>o} + I{a:=o}- (35) 

0 — 

Taking (34) and (35) into account, we come back to (33) and send 0 to 1. Fatou’s lemma 
used to the ds-integral on the left-hand side of (33) and dominated convergence used to the 
rest integrals give 

+ r x>o}{{SYs)+r-^\szs\^ds 

^ JtAr 

< {{SYXr+P r X> 0 }im)+r~\FAs,Ys,Z,) - FXX,Zi) + Xdys)^ +l 2 \SZA)ds 

J tAr 

-p [ X>0}{{SYs)+r~^5ZsdWs. (36) 

J tAr 
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Moreover, (A.2)(ii) implies 


IiSY.>o}{Fi{s,Y,,Zs) - Fi{sX,Zi)) <hsY^>o}{l3i{SYs)+ +ji\6Zs\). 

We then use this inequality to (36). To eliminate the local martingale, we replace r by its 
localization sequence {T„}„gp}+. By the same way of estimation as in Lemma 7 (a priori 
estimate (ii)), we obtain 


{{6Yt^rJ+r <cE[{{6YrJ+r\Ft], 

where c is a constant only depending on T,/3i,/ 32 , 71 , 72 ,^- Since Y,Y' G and P-a.s. 

< Cl dominated convergence yields P-a.s. 1* < Y/. Finally, by the continuity of Y and 
Y' we conclude that P-a.s. Y. < Y!. 


□ 

As a byproduct, we obtain the following existence and uniqueness result. 

Corollary 10 (Uniqueness) Let {A. 2 ) hold for {F, C ■ Then there exists a unique solution 
in SP X MP. 

Proof. (A.2) implies (A.l). Hence existence result holds. The uniqueness is immediate 
from Theorem 9 (comparison theorem). 


□ 

It turns out that a stability result also holds given the convexity condition. We denote 
satisfying (A.2) by {F,Fi,F 2 ,f). We set := N+ U {0}. 

Proposition 11 (Stability) Let p > 1 . Let (A”, F", ^")„gf^o satisfy (A. 2 ) associated 

with (a",/ 3 i,/32,7i, 72) (/J,/), and (Y",Z^) be their unique solutions in SP X A 4 p, respec¬ 
tively. If fC ~ C —^0 and Jq \F^ — F^\{s,Yf^, Zg)ds — 5>0 in as n goes to -|-cxd, then 
(T"’,Z") converges to {Y, Z) in Sp x Mp. 

Proof. We prove the stability result in the spirit of Theorem 9 (comparison theorem). For 
any 6 G (0,1), define 


SFf := F°(t, YC, Z°) - F^{t, Y^, Z?), 

SgY^ := yO - 6»y”, 

SY^ := Y° - r”. 


and 6eZ'^,5Z^, etc. analogously. We observe the 0-difFerence of the generators. Likewise, 
(A.2) (hi) implies that 


F°{t,YC,Z°)-9F^{t,YC,Zf) 

= SFff + (F"(t, Y°, yO) - 0F”(t, r*", Zf)) 


< SFC + {Ffit, Y°, Z°) - OFfit, YC, Zf )) + (1 - e)F^{L 


SgYC SgZl 


i-e’ 1-1 


We first prove convergence of T" and later use it to show that Z" also converges. 
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(i). By exactly the same arguments as in Theorem 9 but keeping along the deduc¬ 
tions, we obtain 


mrry 


pip - 1 ) r 




< +p i{SYp>,}m:)+Y-W5F:\ + {p,+p^y5y:)+ + {^^+^2)\5z:\)ds 

-pl^l[SY^^>o}iiSYY)+Y-^SZ:dWs. (37) 

By the same way of estimation as in Lemma 7 (a priori estimate (ii)), we obtain 

E[(((5F")+)*)n <c(E[(((5r)+)n +E[||57^."||^]), 

where c is a constant only depending on T,/3i,/32,7i, 72 jP- Interchanging and F” and 
analogous deductions then yield 

E[(((-5F-)+)*)n <c(E[((-<5r)+)n +E[||5F."||^]). 

Hence a combination of the two inequalities implies the convergence of F". 

(ii). To prove the convergence of F", we combine the arguments in Lemma 6 (a priori 
estimate (i)) and Theorem 9. To this end, we introduce the function v defined in Section 2 
associated with a function of class X to be determined later. By Ito-Krylov formula, 

v{6eYY) = v{SsC) + r v'{5eY:){F\sX.Zl) - eF-{s,YY, Z:))ds 

Jo 

^"i^oYr)\d0Z:\^ds- vY6eYY)SeZyWs. (38) 

Note that (A.2)(ii)(iii) and v'(SffYY) = sgn((56iF")|u'((5eFg”)| give 

v'(S,YY)(fYs,YY FO) - 0F"(s,YY, F-)) 

< |F(<i,f,")||5f;| 

+ lv'(S,YY)lsgn(S,YY)(FryYYZ°)-0Fr(s,YY,Z:)) 

+ |F(<5,F,")|(/32|<5.F-| +72|5eF;| + Y SoZy). (39) 


We associate v with 


f( I I ) 

so as to eliminate the quadratic term. Note that 
lim v(x) = -IxP, 


lim v'{x) = X. 


0^1 


(40) 


With (39), (40) and (A.2)(ii), we come back to (38) and send 0 to 1. This gives 


\szyds < -i< 5 rp + / + (iai + imsY^i + y+^,)\sz:\)ds 


SYrSZ'^dWs. 
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Now we use the same way of estimation as in Lemma 6 to obtain 


<cE[((,5y”)*)P+||5F."||P], 

where c is a constant only depending on T,/3i,/ 32 , 71 , 72 )P- The convergence of Z” is then 
immediate from (i). 


E 


\5Z:\^ds 


□ 

Remark. So far we have obtained the existence and uniqueness of a LP(p > 1) solution. 
The solvability for p = 1 is not included due to the failure of Lemma 7 (a priori estimate 
(ii)). One may overcome this difficulty by imposing additional structure conditions as in 
Briand et al [3], [6]. To save pages the analysis of solutions is hence omitted. 

6 Applications to Quadratic PDEs 

In this section, we give an application of our results to quadratic PDEs. More precisely, 
we prove the probablistic representation for the nonlinear Feymann-Kac formula associated 
with the BSDEs in our study. Let us consider the following semilinear PDE 

dtu{t, x) + £u(t, x) + F(t, X, u{t, x), cr^ Va;u(t, x)) = 0, 

u{T,-)=g, (41) 

where C is the infinitesimal generator of the solution to the Markovian SDE 

Xt=XQ+ f b{s,Xs)ds+ f a{s,Xs)dBs, (42) 

Jto Jto 

for any (to,xo) G [0,T] x K”, t G [to,T]. Denote a solution to the BSDE 

Yt = g{X*^°’^°) + F{s,Xl«’-\Y,,Z,)ds-Z,dW,, t G[h,Tl (43) 

by (y*o>^ci, 2'*°’“°) or (T, Z) when there is no ambiguity. The probablistic representation for 
nonlinear Feymann-Kac formula consists of proving that, in Markovian setting, u(t, x) := 
is a solution at least in the viscosity sense to (41) when the source of nonlinearity F 
is quadratic in Xxuit, x) and g is an unbounded function. To put it more precisely, let us 
introduce the FBSDEs. 

The Forward Markovian SDEs. Let b : [0,T] x R” —>• R", a : [0,T] x R'^ ^ R”^"^ 
be continuous functions and assume there exists /3 > 0 such that P-a.s. for any t G [0,T], 
|&(t, 0)1 -I- |(T(t, 0)1 < {3 and b{t,x),a{t,x) are Lipschitz-continuous in x, i.e., P-a.s. for any 
t G [0, T], X, x' G R", 


\b{t, x) — b{t, x')\ + \o'{t, x) — a{t, x')\ < /3|x — x'\. 

Then for any (to,Xo) G [0,T] x R”, (42) has a unique solution X*°’^° in for any p > 1. 

The Markovian BSDE. We continue with the setting of the forward equations above. 
Set q > 1. Let Fi,F 2 : [0,T] x R" x R x R'^ —> R, p : R" —>• R be continuous functions, 
(fi : R+ —>■ R+ a continuous nondecreasing function with p(0) = 0 and /(| • |) G I, and 
assume moreover F = Fi + F 2 such that 
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(i) Fi (t, X, y, z) is monotonic in y and Lipschitz-continuous in z, and F2 (t, x, y, z) is mono¬ 
tonic at 1 / = 0 and of linear-quadratic growth in z, i.e., for any {t,x) £ [0,T] x R", 
y,y' e R,z,z' € 

sgn(y - y'){Fi{t, X, y, z) - Fi{t, x,y',z)) < P\y - y'\, 

\Piit,x,y,z) - Fi{t,x,y,z')\ < / 3 |z - z'j, 

sgn(y)F2(t,x, y, z) < / 3 |y| -f / 3 |z| -f /(|y|)|2;p; 

(ii) (y, z) I—S> F2{t, x, y, z) is convex ; 

(iii) for any {t, x, y, z) € [0, T] x R” x R x R'^, 

\F{t,x,y,z)\ < P{1 + |a;|« -f 2|z|) F ip{\y\) F f {\y\)\z\^, 

|y(a;)| < /3(l -f |a;|«). 

Since for any p > 1, the above structure conditions on F and g allow one 

to use Corollary 10 to construct a unique solution in x of (43) for 

any p > 1. Moreover, by standard arguments, is deterministic for any (to,xo) £ 

[0, T] X R". Hence u{t, x) defined as is a deterministic function. With this fact we now 
turn to the main result of this section: u is a viscosity solution of (41). Before our proof 
let us recall the definition of a viscosity solution. 

Viscosity Solution. A continuous function u : [0,T] x R" —)• R is called a viscosity 
subsolution (respectively supersolution) to (41) if u{T,x) < g{x) (respectively u{T,x) > 
g{x)) and for any smooth function (j) such that u — (j) reaches the local maximum (respectively 
local minimum) at (to, 2 : 0)1 we have 

dt(l){to,Xo) FC<j){to,xo) F F{to,xo,u{to,xo),cr^yx(t){to,Xo)) > 0 (respectively < 0). 

A function u is called a viscosity solution to (41) if it is both a viscosity subsolution and 
supersolution. 

Proposition 12 Given the above assumptions, u(t,x) is continuous with 

\u{t,x)\ < c(l -f |x|«), 

where c is a constant. Moreover, u is a viscosity solution to (41). 

Proof. Due to the Lipschitz-continuity of b and a, is continuous in (t,x), e.g., in 

mean square sense. The continuity of u is then an immediate consequence of Theorem 11 
(stability). The proof relies on standard arguments and hence is omitted. By Lemma 7 (a 
priori estimate (ii)), we prove that u satisfies the above polynomial growth. It thus remains 
to prove that u is a viscosity solution to (41). 

Let (/) be a smooth function such that u — (j) reaches local maximum at (to, xo). Without 
loss of generality we assume that the local maximum is global and u(to, xq) = (/)(to, xq). We 
aim at proving 

dt(l)(to,xo) F C(l)(to, Xo) F F(to, xo, u(to, xq), a^Vx(l)(to,xo)) > 0. 

From (43) we obtain 

Yt=Yt,- f F(s,Xl°^^°,Ys,Zs)dsF f Z^dWs- 

Jto dto 
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By Ito’s formula, 


=0(fo,a:o)+ f {dsct> + Cct>}{s,Xlo’^°)ds+ f 

Jto jtQ 

Now we take any t G [tQ,T]. Note that the existence of a unique solution of (42) and (43) 
implies by Markov property that Yt = Hence, (j){t,Xl°'^°) > u{t,Xl°'^°) = Yf 

By touching property, on the set Xl°"^°) = k^t} we have 

dtct>{t,Xl°’^°) + £<j){t,Xl°’^°) + F{t,Xl°’^°,Yt,Zt) >0 P-a.s., 
a'^Va:ct>{t,Xl°’^^)- Zt = 0 P-a.s. 

Now we set t = to. We have (j){to,Xl°'^°) = (j}{to,xo) = u{to,xo) = Ytg. Moreover, the 
above equality implies Ztg = cr^Va;(/)(to, a^o)- Plugging the two equalities into the above 
inequality gives 

dt(l){to,xo) -I- C(l)(to,xo) + F{to,xo,u{to,xo),cr^Xx(l){to,xo)) > 0. 

Hence u is a viscosity subsolution of (41). u being a viscosity supersolution and thus a 
viscosity solution can be proved analogously. 


□ 
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